The effect of irradiation-induced disorder on the conductivity and critical 
temperature of the organic superconductor k-(BEDT-TTF) 2 Cu(SCN) 2 
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We have introduced defects into clean samples of the organic superconductor k-(BEDT-TTF)2- 
Cu(SCN)2 in order to determine their effect on the temperature dependence of the conductivity and 
the critical temperature T c . We find a violation of Matthiessen's rule that can be explained by a 
model of the conductivity involving a defect-assisted interlayer channel which acts in parallel with 
the band-like conductivity. We observe an unusual dependence of T c on residual resistivity which is 
not consistent with the generalised Abrikosov-Gor'kov theory for an order parameter with a single 
component, providing an important constraint on models of the superconductivity in this material. 

PACS numbers: 74.70.Kn, 74.62.Dh, 74.25.Fy 



Quasi-two-dimensional organic conductors based on 
the donor molecule bisethylenedithiotetrathiafulvalcnc 
(BEDT-TTF or ET, see inset to Fig. 1) have attracted 
sustained interest Q . These materials consist of layers of 
ET separated by inorganic anions and can be prepared 
with exceptional purity, so that magnetic oscillations in 
the resistivity can be observed at relatively low fields . 
Many experiments have been conducted to elucidate the 
symmetry of the superconducting order parameter but 
despite growing evidence of unconventional superconduc- 
tivity [J, the matter has remained unresolved Q. 

The normal state of these materials is unusual with 
the temperature dependence of the electrical transport 
and the optical conductivity deviating significantly from 
what would be expected for a conventional metal , 
showing a crossover from insulating-like conductivity at 
high temperature to metallic-like behaviour at lower tem- 
perature (see Fig. 1). This effect has been successfully 
described within the framework of dynamical mean-field 
theory (DMFT) as a crossover from an incoherent "bad- 
metal" state at high temperatures to a coherent Fermi 
liquid below ^30 K [jj |g. The usual effect of disorder 
is to change only the temperature independent compo- 
nent of the resistivity; this is known as Matthiessen's 
rule. However, several strongly correlated materials vio- 
late Matthiessen's rule including various cuprate |7j and 
organic || superconductors. Violations of Matthiessen's 
rule have previously been taken as evidence for non- Fermi 
liquid behaviour. In particular, Strack et al. || have ar- 
gued that the violations of Matthiessen's rule in the salt 
K-(ET) 2 Cu[N(CN) 2 ]Br indicate that the description of 
the conductivity in terms of a crossover from an incoher- 
ent metal to a Fermi liquid is incorrect. 

The behaviour of the superconducting transition tem- 
perature T c as a function of non-magnetic disorder 
lJ E3- El also provides a crucial test of the symmetry 



of the order parameter. Anderson's theorem ^lj states 
that for s-wave pairing non-magnetic impurities do not 
affect T c . For magnetic impurities T c is strongly re- 
duced for singlet states in a manner described by the 
Abrikosov-Gor'kov (AG) formula For non- s-wave 

(including extended-s) order parameters, scattering due 
to non-magnetic impurities reduces TV in a manner again 
described using the AG formula jj, |fj, [lfj . 

In order to investigate these effects experimentally for 
the ET superconductors, we have selected k-(ET)2Cu- 
(SCN)2, which has one of the highest transition temper- 
atures (T c « 10 K) in this family and whose electronic 
properties have been the subject of detailed study 0]. 
It is possible to introduce disorder into these samples 
by adjusting the cooling rate 0, ^| , an effect which is 
probably due to the freezing of terminal ethylene group 
disorder which occurs at around 70-80 K [J, [lfj . How- 
ever, we find that we have more control in our experiment 
by introducing disorder via irradiation by either x-rays 
or protons. X-rays and protons give similar effects, but 
we have been able to perform more detailed experiments 
using the former and hence we concentrate mainly on the 
x-ray damaged samples in our discussion. 

In this Letter we report the violation of Matthiessen's 
rule in k-(ET)2Cu(SCN)2- However, we find that a sim- 
ple theory which includes the effect of interlayer scatter- 
ing from impurities is able to explain this effect. Further 
we have confirmed quantitative predictions of this model 
relating the low temperature resistivity to the high tem- 
perature conductivity. Thus we find that the violation 
of Matthiessen's rule is consistent with the DMFT de- 
scription of the crossover in transport behaviour |fj] . We 
also observe a disorder-induced suppression of T c which 
is consistent with a non-s-wave gap. 

Our measurements have been performed on single crys- 
tals of k-(ET)2Cu(SCN)2 grown electrochemically. Gold 
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wire of thickness 12.5 /Ltm was attached to the sam- 
ples (typical dimensions 0.3x0.3x0.1 mm 3 ) with graphite 
paste in the four-probe configuration. Samples were 
cooled from 120 K to 10 K at a rate of 20 K/hr, below the 
slowest rate used in Ref. [lj, in order to avoid introducing 
disorder from fast-cooling. Defects in our samples were 
created at room temperature by filtered Cu K a radiation 
(E = 8 kcV, A = 1.54 A) from a Cu x-ray rotating anode 
(typically 55 kV, 50 mA) at a flux of « 2.5 x 10 s photons 
s — 1 . The computer program RADDOSE [lj| allowed this 
to be converted into a dose rate of « 105 Gy s _1 . Using 
tabulated absorption coefficients and assuming that the 
mass absorption coefficient of k-(ET)2Cu(SCN)2 can be 
calculated as the sum of the mass absorption coefficients 
of the constituent atoms 0]> we estimate the x-ray at- 
tenuation length to be 90 /xm, approximately the thick- 
ness of the samples. In order to attain uniform damage 
we irradiated both sides of the sample. X-ray doses up 
to 630 MGy were used. Proton irradiation experiments 
took place at Surrey Ion Beam Centre. Protons were ac- 
celerated to 4 MeV and could be implanted to a mean 
depth of 150 /xm, providing an approximately uniform 
damage profile for samples <100 /im. For each incre- 
mental radiation dose, we make a measurement of the 
transport properties; the contact configurations stay the 
same throughout the experiment. For each type of irra- 
diation, the resistivity was found to be reproducible over 
multiple thermal cycles, so that we can be confident that 
the observed changes are due to irradiation and not to 
thermal cycling. 

It is immediately apparent from Fig. 1 that the effect 
of increasing irradiation dose is to decrease the resistivity 
over most of the temperature range, and in particular to 
reduce the magnitude of the broad peak centred around 
T p ^90 K. Well into the Fermi- liquid regime however, the 
behaviour recovers a traditional metallic character; below 
^46 K, the resistivity increases with increasing irradia- 
tion dose (see inset to Fig. 1 and Fig. 2). This violation 
of Matthiessen's rule is not predicted by DMFT alone 
for the low defect densities produced in our experiments 
and therefore appears at first sight to be at odds with 
the DMFT description of electronic transport in layered 
organic charge transfer salts. 

The low-temperature region of Fig. 1 is presented in de- 
tail in Fig. 2, more clearly showing the superconducting 
transition for the same range of x-ray irradiation doses. 
For each trace, the resistivity follows a Fermi-liquid like 
T 2 dependence at temperatures above the superconduct- 
ing transition and a sharp drop at the onset of super- 
conductivity. The residual resistivity po is extracted by 
fitting data to p — po + AT 2 in the temperature range 
from just above T c to 20 K. The upper inset to Fig. 2 
presents the change in residual resistivity Apo with re- 
spect to an undamaged sample as a function of dosage, 
showing an approximately linear dependence. 

The observation that increasing defect density in- 




t i;k; 



FIG. 1: Interlayer resistivity of an x-ray damaged sample as 
a function of temperature for a range of irradiation doses, 
showing the violation of Matthiessen's rule. Data from the 
proton irradiated samples are qualitatively the same. Insets 
show the ET molecule and a magnified view of the region near 
Tcroas (defined in the text). 
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FIG. 2: The low temperature resistivity of k-(ET) 2 Cu(SCN) 2 
increases with x-ray dose in an approximately linear fashion, 
as shown in the upper inset. Data from proton-irradiated 
samples show a similar pattern as a function of disorder. The 
lower inset shows three methods of defining the critical tem- 
perature: a lower T c i, an upper T C 2 and T c , the maximum of 
the first derivative, which is the definition used in this paper. 
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creases the interlayer conductivity over a wide temper- 
ature range leads us to hypothesize that defects affect 
the conductivity in two ways: (i) the resistivity associ- 
ated with band-like transport due to the overlap of the 
molecular orbitals increases linearly with defect density 
as prescribed by Matthiessen's rule; (ii) there is a parallel 
defect-assisted interlayer channel, whose conductivity is 
proportional to the defect density 0]. This model sug- 
gests that the interlayer conductivity <r(x, T) depends on 
a dimcnsionless quantity x, which is proportional to the 
defect density, and temperature T as 



a{x,T) 



1 



Po(0) + •^Pimp /^intrinsic 

(T) 



+ xa±, 



(1) 



where po(0) is the residual resistivity of the undam- 
aged sample, xpi mp is the contribution to the resis- 
tivity from defect scattering in the transport due to 
molecular orbital overlap, Pmtrinsic(^) is the intrinsic 
temperature-dependent scattering contribution (due to 
electron-electron interactions, phonon scattering, etc.) 
and xa± is the defect-assisted interlayer conductivity. 

The applicability of this model can be demonstrated 
by examining low and high temperature limits. At low 
temperature pintrinsic(T) is small so that 



*£/?imp i ^intrinsic 



(2) 



and hence the resistivity increases linearly with defect 
density, x. At high temperature Pintrinsic(T) 3> Po(0) + 
a^Pimp and Eqn Q becomes 



p(x,T) w (xa± + p. m l liasic (T)) 



(3) 



hence the conductivity p(x, T) 1 linearly increases with 
x, as observed. The inset to Fig. 1 shows that there is 
a temperature, T cross , at which the resistivity is inde- 
pendent of defect density. In our model, this can be 
found by evaluating da(x, T)/dx — 0. We find that 
a T-independent crossing point occurs when the condi- 
tions (i) Xp imp < p (0) + ^intrinsic 

and (ii) p(0,T cross ) = 
(Pimp/ojJ 1 / 2 are satisfied. For our x-ray damaged sam- 
ples, this yields an estimate of pi mp /crx = 3x 10 3 £l 2 cm 2 . 
We can obtain an independent estimate of this pa- 
rameter by plotting Aa p = a(x,T p ) — er(0, T p ) against 
Apo = Po(x) — Po(0), as shown in Fig. 3. The for- 
mer quantity can be evaluated in the high temperature 
limit (Eqn. El to be Aa p « xaj_, while the latter quan- 
tity can be evaluated in the low-temperature limit to be 
Apo ~ 2:pimp- The approximate straight-line dependence 
(which breaks down when Apo and Aoo are large) ob- 
served in Fig. 3 shows that both quantities are param- 
eterised by x in the same way, and the gradient of the 
line yields pim P /cr± = 1.5 x 10 3 f2 2 cm 2 , in order of mag- 
nitude agreement with our previous estimate [lfij . We 
note that for proton irradiation, both the gradient of the 
line in Fig. 3 and the temperature T cross are increased, 
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FIG. 3: The dependence of the change in conductivity at the 
peak, Act, on the change in residual resistivity, Apo, which is 
approximately linear, except for at large Apo. Due to the dif- 
ferent nature of the defects, the slope of proton (two different 
samples shown with open circles and squares) and x-ray (one 
sample indicated by filled circles) data differ, though the ef- 
fect remains qualitatively the same. Inset: the position T p of 
the resistivity peak shifts slightly to higher temperature with 
increasing po(x). 



demonstrating that the nature of the irradiation affects 
the ratio pim P /a±, suggesting that x-rays and protons 
produce different types of damage (protons producing a 
more effective interplane transport channel than x-rays). 

The value of T c , corresponding to the maximum of 
dp/ 'AT in Fig. 2, is plotted as a function of po in Fig. 4. 
We find that T c falls with defect density (in agreement 
with measurements of defects induced by cooling-rate dis- 
order |l3|h but the dependence exhibits a sharp change in 
gradient when po reaches a threshold value Pq ~ 2 17cm. 
However, even at the highest defect densities studied, 
the samples exhibit a superconducting ground state. We 
also find that the width AT = T c2 - T cl (see Fig. 2) of 
the superconducting transition decreases with increasing 
defect density, as shown in the inset to Fig. 4. This is 
consistent with T c exhibiting a change of gradient with 
damage; at high damage, local variations in damage have 
a smaller effect on the broadening of the transition be- 
cause dT c /dpo is lower. 

The theory of AG H,E3 

for non-magnetic defects in a 
non-s-wave superconductor implies that the suppression 
of T c follows the AG formula given by 



J-cO \2 



' 2 + 47rfc B T c r 



(4) 



where tp is a digamma function and r is the scattering 
time. In the low-defect density limit (po < Pg), this 
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radiation-induced disorder in k-(ET)2Cu(SCN)2- We 
find that a dramatic departure from Matthiessen's rule 
can be straightforwardly explained in terms of defect as- 
sisted interlayer tunnelling. Although T c initially follows 
a dependence on po consistent with pair-breaking scatter- 
ing, the superconducting state proves to be robust into 
the dirty limit. This unusual dependence of T c on po 
is not consistent with an order parameter with a single 
component, providing an important constraint on models 
of the superconductivity in this material. 
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FIG. 4: The suppression of critical temperature T c compared 
to the residual resistivity po(x). The suppression is approxi- 
mately linear for low damage po(x) < po = 2 f!cm, but there 
is a change in gradient for po > Pa- The inset shows the de- 
pendence of the width AT C of the superconducting transition 
on po(x). Symbols correspond to those in Fig. 3 and data 
for an additional undamaged sample are shown by a fillled 
diamond. 

yields T c q — T c ~ irfi/^k-QT. This linear sector should 
have a slope consistent with interlayer transport theory 
0, i-e. dT c /dpo = -e 2 m*d_L^/(4fc B ?i 3 ), where m* is the 
effective mass, d± is the interlayer spacing and t± is the 
interlayer transfer integral. Using the value of m* from 
transport measurements, our measured slope of the sup- 
pression of T c for po < Po yields (the line in Fig. 4), a 
value of t± = 0.03 ± 0.01 meV, in good agreement with 
the value from angle-dependent magnetotransport [20| . 
Thus the suppression of T c for po < Po is consistent with 
a non-s-wave gap. However, the departure from the AG 
formula for po > Pq casts doubt on this interpretation. 

We note that our value of pg corresponds to a ratio 
of the in-plane coherence length to the mean free path 
£/Z « 0.2. If all scattering events contributing to I were 
pair-breaking, it would be expected that superconduc- 
tivity should be suppressed by this degree of scattering. 
Crucially, the observed change of gradient is not expected 
for a model involving an order parameter with a single 
component within the AG theory. Even a generalised 
AG equation describing multicomponent order param- 
eters | 2l| does not quantitatively agree with our data, 
though this approach assumes that disorder does not af- 
fect the symmetry of the order parameter; it is probably 
necessary to examine specific candidate pairing interac- 
tions in detail. An explanation for this behaviour might 
involve a mixed order parameter with both s-wave and 
unconventional components or interband scattering. 

In conclusion, we have investigated the effect of 
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